Time-Resolved Spin Torque Switching and Enhanced Damping in Py/Cu/Py 

Spin- Valve Nanopillars 



in 
o 
o 

> 
O 



N. C. Emley, I. N. Krivorotov, A. G. F. Garcia, O. Ozatay, J. C. Sankey, D. C. Ralph, and R. A. Buhrman 

Cornell University, Ithaca, NY 14853-2501 
(Dated: February 2, 2008) 

We report time-resolved measurements of current-induced reversal of a free magnetic layer in 
Py/Cu/Py elliptical nanopillars at temperatures T = 4.2 K to 160 K. Comparison of the data to 
Landau-Lifshitz- Gilbert macrospin simulations of the free layer switching yields numerical values 
for the spin torque and the Gilbert damping parameters as functions of T. The damping is strongly 
T-dependent, which we attribute to the antiferromagnetic pinning behavior of a thin permalloy 
oxide layer around the perimeter of the free layer. This adventitious antiferromagnetic pinning layer 
can have a major impact on spin torque phenomena. 

PACS numbers: 85.75.-d, 75.75.+a, 81.65.Mq 
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Experiments Q, H, H, S have shown that a spin- 
polarized current passed through a nanomagnet can ex- 
cite a dynamic response as the result of a spin torque 
applied by the conduction electrons |^ . The potential 
for technological impact of this spin transfer (ST) effect 
has inspired research in DC current-induced microwave 
oscillations 0,0 and hysteretic switching 0,001 in cur- 
rent perpendicular to the plane (GPP) nanopillars and 
nanoconstrictions. Typically, ST switching data is ob- 
tained through the use of slow current ramp rates ('--^l 
mA/s), but fast pulses (^^10^° mA/s) access the regime 
where thermal activation of the moment over a current- 
dependent barrier ^,^9,1 does not play a major role in the 
switching process. This spin torque-driven regime [loj 
is advantageous for the quantitative examination of the 
spin torque parameters due to the computational acces- 
sibility of numerically integrating the Landau-Lifshitz- 
Gilbert (LLG) equation for short durations. 

Here we report time- resolved measurements of the spin 
torque-driven switching event in Gu 100/Py 20/Gu 6/Py 
2/Gu 2/Pt 30 (in nm, Py = NigiFeig) GPP spin-valve 
nanopillar structures at bath temperatures T = 4.2 K to 
160 K. We compare our experimental results with LLG 
simulations in the macrospin approximation and find 
good agreement between simulation and measurement. 
This both confirms the applicability of the macrospin ap- 
proximation in the spin torque-driven regime and facili- 
tates the quantitative determination of T-dependent spin 
torque and magnetic damping parameters. At higher T 
we find that the strength of the spin torque exerted per 
unit current is in reasonable numerical accord with re- 
cent model calculations, and that the damping parameter 
aQ for the nanomagnet excitations is both anomalously 
high, as suggested by previous pulsed current measure- 
ments rnj, and T-dependent. The strong T variation of 
ao, in conjunction with anomalous behavior of the nano- 
magnet switching fields Hs,i{T) in some devices, points 
to the presence of an adventitious antiferromagnetic ox- 
ide layer around the perimeter of the nanomagnet that 
has a major effect on the nanomagnet dynamics driven 
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FIG. 1: (a) Slow ramp rate spin torque switching of a nano- 
magnet as measured by the GMR effect for sample 1, a 



60x190 nm ellipse, at T = 40 K and 



385 Oe, which 



opposes the dipole field so that -ffapp -I- -ffdip ~ 0. Arrows 
indicate the scan direction, (b) Pulsed spin torque P to AP 
switching measured for the same sample at / = 1.07 mA (□) 
and 2.13 mA (■). The data (symbols) have been normal- 
ized to Mx = ±1 for simple comparison with the simulated 
macrospin switching (lines). Pulse shape distortions are due 
to the setup (see ref. [13]). 



by a spin torque. 

The nanopillar devices employed in this study were 
fabricated using a process described elsewhere A 
slow ramp rate ST scan is shown in Fig. 1(a) for sample 
1, a 60 X 190 nm ellipse. In Fig. 1(b) we show parallel (P) 
to anti-parallel (AP) switching events for sample 1, aver- 
aged over 10,000 switches, taken at pulsed current ampli- 
tudes / = 1.07 mA and 2.13 mA at T = 40 K as open and 
solid squares, respectively. The measured signal is a time- 
resolved voltage drop \I ■ A_R| from the giant magnetore- 
sistance (GMR) of the sample as the free layer switches 
from P to AP orientation, where Ai? = i?3;(AP)-i?2:(P) 
and Rx is the 4-point device resistance. The data have 
been normalized to = -\- / — \ (minimum resistance / 
maximum resistance) for simple comparison with simu- 
lated switching events, described below, which are shown 
as solid and dashed hues. The abrupt ('^200 ps) jump 
from Mx = -1 to 1 at time = is not a switching event but 
is simply the rising edge of the current pulse. The more 
gradual transition between P [M^ — 1) and AP (return 
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to Mx = -1) is the envelope coming from averaging over 
thousands of individual switching events, each of which 
follows a trajectory determined by initial conditions that 
are randomized by the stochastic thermal fluctuations of 
the free layer. We define the switching time tswitch as 
the time elapsed between 50% of the signal rise and 50% 
of the signal drop as indicated in Fig. 1(b) 

To obtain a quantitative understanding of the ST 
switching, we have simulated the nanomagnet dynam- 
ics by numerical integration of the LLG equation in 
the macrospin approximation with the inclusion of a 
Slonczewski-type spin torque term. 

^ = 7[m X [HcS + -f^Lang(T')) - a{e)m x (to x (i?off 

- eM..(TO?afJ.d)sin. - ^ ^ ^ 

Here 7 is the gyromagnetic ratio, to is the unit direc- 
tional vector of the free layer macrospin, p is the spin 
polarization axis, is the in-plane angle between them, 
g{6) is the spin torque function, Ms{T) is the free layer 
magnetization, as measured separately for a continuous 
2 nm Py film in a Cu/Py/Cu trilayer that was exposed 
to the same heat treatments as the nanopillars, d is the 
nanomagnet thickness, area = is its lateral area with 
dimensions a and h that are estimated by OOMMF mi- 
cromagnetic simulations (see below), and HeS is the 
sum of external 7?oxt, in-plane anisotropy i?Kj and out- 
of-plane anisotropy H± fields. T^oxt is the sum of the 
magnetostatic dipole field from the fixed layer i/dip fmd 
the applied field i?app from the electromagnet, which is 
adjusted to compensate for _ffdip so 7?cxt ~ 0. 

The initial conditions of the simulation were set by 
0^ = + 0^is + ^randCT), whcrc - 0° for P to 
AP and 180° for AP to P switching and 6'mis repre- 
sents any systematic angular misalignment between free 
and fixed layer moments due to the setup and was gen- 
erally set to 0. The random angle 0rand(r) is treated 
as a Gaussian with a standard deviation y^kBT/2EQ{T) 
where Eo{T) = Eq{A.2 K)[Ms(T)/M,(4.2 K)]^ is the uni- 
axial anisotropy energy. Ea{A.2 K), a, and b are esti- 
mated from T = 0, 2D OOMMF simulations of Py el- 
liptical disks having Hk and AR values similar to those 
measured at 4.2 K. The lateral area is estimated to 12% 
uncertainty with this method, nearly a factor of 2 bet- 
ter than the inherent shape variation among otherwise 
identical elliptical patterns due to lithographic fluctua- 
tions. Ohmic heating effects during the current pulse are 
taken into account by locally raising the temperature of 
the device to T' = ^T^ + 10.23(K/mV)2(i?^(r) • 7)2 0. 
A Langevin field Hha.ng{T') accounts for thermal fluc- 
tuations during the dynamic trajectory, fluctuating 
randomly in 3-dimcnsions with a standard deviation 
^y2aokBT'iiohMJ'r){area ■ d)At where At = 1 ps is 
the time step [s], [lOl ll5j . Gilbert damping is assigned an 
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FIG. 2: Measured (symbols) and simulated (lines) 1/t switch 
versus / for sample 2, an 80x180 nm ellipse, at T = 160 K 
(□), 40 K (•), and 4.2 K (a) for (a) AP to P and (b) P to 
AP switching. Simulations to 1/t switch ~ yield estimates of 
the intercepts /^^(T). (a) inset: I^g{T)/M^{T). 

angular dependence a(6') — aoll — i/sin^ 6/{l — i'^ cos^ 9)], 
where v = 0.33 for Py/Cu/Py nanopillars but the 
addition of this angle-dependent damping term had only 
a small effect on the simulation results. 

The spin torque function is approximated by g{9) — 
A sin 6'/(l -|- B cos 9) where A and B are phenomenologi- 
cal parameters |lllll7lll8l |. In our simulations we use ag, 
A, and B as T-dependent fitting parameters to match 
the simulated with the measured values of 1/ 1 switch ver- 
sus / for each T, where we allow to be different for 
the two switching directions. In Fig. 1(b) we plot the 
average of 2000 simulated P to AP switching events at 
T = 40 K alongside the normalized data for sample 1 
with the best fit simulation yielding A = 0.5, B = 0.11, 
and oiQ = 0.048. Since the current step in the simu- 
lation turns on instantaneously, an average pulse half 
rise time of 112 ps, measured from data such as those 
in Fig. 1(b), has been added to all simulated t switch- 
We plot measured l/tswUch versus / for AP to P and 
P to AP switching at T = 160 K, 40 K, and 4.2 K for 
sample 2, an 80x180 nm ellipse, together with best fit 
simulations, all of which are averages over 2000 events, 
in Fig. 2(a) and 2(b), respectively. Simulations out to 
long switching times {l/tswitck < 0.1 ns~^) allow for 
good estimates of the l/tswitck ^ intercepts /^(T), 
which are the critical currents (+ = P to AP) defin- 
ing the onset of spin torque-driven switching. These 
should depend on the spin torque and damping pa- 
rameters as I^{T) oc aoM'^{T) 10]. A striking result 
from these measurements is the strong T-dependence of 
lf^[T)/M'^{T) (Fig. 2(a) inset), which varies by more 
than 60% over the entire T range, where the upturns 
at low T indicate a strong dependence of damping, spin 
torque, or both. 

In Fig. 3 the best fit values for ao. A, and B (assum- 
ing 6'mis = 0°) are plotted as functions of T for sample 2. 
Uncertainties in the fit parameters, Aao = 0.0035, = 
0.025, and AB = 0.045, are found through an exploration 
of parameter space about the best fit values. Account- 
ing for these T-dependences, the theoretical prediction 
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FIG. 3: Best fit parameters (a) damping ao (for P to AP (o) 
and AP to P (•) switcliing) and (b) spin torque parameters A 
(a) and B (a) as functions of T from matciiing simulated with 
measured values of l/tswHch versus / for each T for sample 2. 



of I%{T) oc a^(T)M2(r)(l ± B{T))/AlT) agrees with 
the measurement to within 10% over the entire range of 
T. All four devices that were extensively studied show 
an amplitude and T-dependence of ao very similar to 
that of Fig. 3(a); a gradual but significant increase with 
decreasing T below 160 K, above which the devices are 
thermally unstable, followed by a stronger increase start- 
ing below 60 K - 40 K where the best fit values of also 
suggest differences between the two switching directions. 
For T < 60 K, the trends in the T-dependence of the spin 
torque parameters A and B vary from sample to sample, 
but for T > 60 K both consistently show a very mild 
dependence on T as illustrated in Fig. 3(b). For the four 
samples studied in detail we found that at 40 K A ranged 
from 0.5 to 0.68 and B varied from 0.11 to 0.35. In gen- 
eral we also found that A would decrease by 10 or 20% in 
going from 40 K to 160 K while B would typically vary 
by 10% or less. These values of A and B and the varia- 
tion with T > 60 K can be compared with the results of 
a two-channel model 0| with which the measured GMR 
parameters, R{T) and Ai?(T), can be used to predict 
the spin torque parameters [20|. This model predicts A 
= 0.52, B = 0.36 at 40 K, with A decreasing to 0.47 at 
160 K and B remaining essentially constant. This is in 
reasonable accord with the data, given the experimental 
uncertainties in nanomagnet size and alignment. 

We attribute the significant T-dependence of ao to the 
presence of a weak antiferromagnetic (AF) layer on the 
sidewalls of the nanopillar. Although no such AF layer 
was deliberately deposited, the exposure of the nanopil- 
lars to air after ion mill definition undoubtedly oxidized 
the sidewalls, thus allowing for AF Py oxide to form and 
weakly exchange bias the ferromagnetic layers. An ex- 
ample of direct evidence for this adventitious exchange 



biasing is shown in Fig. 4, where switching fields Hsi 
and Hs2, defined in the inset, from 20 field scans at 
each T, are plotted from 4.2 K to 160 K for sample 3, 
an 80x180 nm ellipse, a previously unmeasured device 
cooled in i?app = 0. Note that Hsi varies more rapidly 
with T than Hs2i which is indicative of an exchange bias 
that strengthens with decreasing T (particularly rapidly 
below 40 K) and promotes AP alignment, i.e. a bias 
set by the dipole field from the fixed layer. Another key 
point illustrated in Fig. 4 is the large variation that de- 
velops in Hsi, and to a lesser extent in Hs2, upon multi- 
ple minor loop scans after the device is cooled to low T. 
Initially, the device switches repeatedly with nearly the 
same switching fields, but after six or seven magnetic re- 
versals the switching fields begin to fluctuate greatly from 
reversal to reversal, indicating stochastic variations in the 
net strength of the oxide pinning field. While the effects 
of the random pinning field are particularly pronounced 
at 4.2 K they are observed up to 160 K, indicating that 
some degree of magnetic ordering within the AF persists 
over this entire T range and also that each reversal of 
the free layer nanomagnet has an irreversible perturb- 
ing effect on the magnetic structure of the AF oxide. It 
is important to note that the slow ramp rate current- 
driven switching events at low T for these devices show 
good reproducibility, with little variation from one sweep 
to another, as should be the case because ST switching 
currents are less sensitive to field variations than are the 
switching fields. The strength of this low T AF exchange 
biasing varies from device to device, with some samples 
showing no random variations in Hs^i ■ We do not believe 
that fluctuations of this sort have affected any previously- 
published conclusions from our group. Nevertheless, the 
fluctuations visible in some samples indicate clearly the 
presence of an AF layer that should influence the prop- 
erties of all nanopillar ST devices. 

Exchange biasing in AF/Py films has been demon- 
strated to dissipate dynamic magnetic energy through 
a two-magnon scattering process arising from local vari- 
ations in the interfacial exchange coupling [21L |23| . 
Over the course of the pulsed / measurements, the free 
layer is switched hundreds of millions of times, which the 
Hs,i data indicate should result in the AF layer being 
on average magnetically ordered but with a finer, more 
randomized local magnetic structure that leads to strong 
damping. The rapid increase in damping, observed over 
the same low T range where both the unidirectional AF 
pinning field and the critical currents (T) also increase 
rapidly, is attributed to an increasing portion of the AF 
oxide layer becoming blocked, thereby simultaneously in- 
creasing the amount of interfacial exchange coupling vari- 
ation seen by the free layer nanomagnet as it moves in 
its dynamic switching trajectory, consistent with the two- 
magnon model. The process of inducing randomization 
in the AF by the nanomagnet reversal itself may also 
lead to enhanced damping. The unidirectional pinning 
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FIG. 4: (a) Switching fields Hs,i of the free layer nanomagnet 
(defined in the inset) measured as a function of T. Hsi (o) 
are AP to P and Hs2 (•) are P to AP switching fields. Pro- 
gression of the random exchange field from the AF perimeter 
layer is observed in subsequent minor loop GMR scans (b) 
through (e). This previously unmeasured device (sample 3, 
an 80 X 180 nm ellipse) was cooled to 4.2 K in -ffapp ~ 0, where- 
upon a single, orientation-setting major loop scan, followed by 
20 GMR minor loop scans were taken. The sample was then 
sequentially warmed back to 160 K in 20 K steps, with 20 
GMR minor loops measured at each T. At 4.2 K, the free 
layer switched consistently at fields Hsi ~ 650 Oe and Hs2 
~ 240 Oe, shown as A and ▲ in (a), respectively, for the first 
six GMR scans. Subsequent scans, however, showed more 
stochastic switching behavior that persisted for the duration 
of the experiment. All 20 Hs,i for each T are shown, although 
some are indistinguishable due to the size of the symbol. The 
solid and dashed lines in (a) are guides to the eye for the 
maximal Hsi and Hs2 values, respectively, as functions of T. 



boundary conditions that should be employed in micro- 
magnetic modeling of the free layer nanomagnet behav- 
ior and affect the dynamical modes of ST-driven preces- 
sion. We are currently investigating whether this AF 
perimeter layer can account, at least in part, for the 
narrower than predicted ST-induced microwave oscilla- 
tor linewidths that have been observed in similar nano- 
magnets at low T |25| . 

In summary, we have performed time-resolved mea- 
surements of the spin torque-driven switching of a Py 
nanomagnet at T = 4.2 K to 160 K. LLG macrospin 
simulations are in close quantitative agreement with the 
ST switching events, yielding values of the parameterized 
spin torque function g{9) = A sin 9 /{I + Bcos9) and the 
damping parameter. We find ao to be high, > 0.03, and 
strongly T-dependent, which we attribute to the AF pin- 
ning behavior of a thin Py oxide layer on the sidewall of 
the nanomagnet. The values of A and B are in fair nu- 
merical agreement with the spin torque calculated from 
the two-channel model using the measured magnetoresis- 
tance values of the nanopillar spin-valve. There is, how- 
ever, considerable device-to-device variation in the spin 
torque asymmetry parameter B, which we tentatively at- 
tribute to the variable nature of the AF perimeter layer. 
The presence of an AF oxide layer can have a major effect 
on the nanomagnet dynamics. Controlling this layer will 
be important in optimizing spin torque-driven behavior. 
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1-0541, and by NSF through the NSEC support of the 
Cornell Center for Nanoscale Systems. Additional sup- 
port was provided by NSF through use of the facilities of 
the Cornell Nanoscale Facility - NNIN and the facilities 
of the CorneU MRSEC. 



field is present over the entire T range, with diminish- 
ing amplitude with increasing T. As the AF grains be- 
come unblocked, an additional damping mechanism be- 
comes possible if these grains can undergo reversal on the 
same time scale that the free layer traces out its dynamic 
switching trajectory, which can result in domain drag or 
the "slow relaxer" dissipation process j^- This effect 
could make a significant contribution to the greater than 
intrinsic damping that persists to higher T. 

Most ST device fabrication processes currently em- 
ployed expose the sides of the free layer nanomagnet to 
some level of an oxidizing ambient at some point, either 
during or after processing and to our knowledge there 
have been no reports of actively protecting the sidewalls 
from oxidation. We suggest that the native AF oxide 
layer that forms can have substantial, previously under- 
appreciated consequences for the ST behavior, leading 
to a substantially enhanced damping parameter which 
directly increases the critical currents for switching. The 
presence of this AF perimeter layer may also alter the 
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